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I Abstract. We investigate Einstein theories of gravity, coupled to a scalar field (f and 

point-like matter, which are characterized by a scalar field-dependent matter coupling 
function e^^'^\ We show that under mild constraints on the form of the potential for 
the scalar field, there are a broad class of Einstein-like gravity models -characterized by 
the asymptotic behavior of H- which allow for a non-Newtonian weak-field limit with the 
gravitational potential behaving for large distances as Inr. The Newtonian term GM/r 
^ , appears only as sub-leading. We point out that this behavior is also shared by gravity 
' models described by f{R) Lagrangians. The relevance of our results for the building 
of infrared modified theories of gravity and for modified Newtonian dynamics is also 
discussed. 

O 

^ '. 1 Introduction 
o 

In the last decade there has been growing interest for gravitational models with large- 
^ ■ dist ance deviations from standard, Einsteinian and Newtonian gravity [H [21 [3l HI [5l [6l [71 
[HI M, HD] • The main reason behind this interest is the hope that the accelerated expansion 
of the universe and rotational curves of the galaxies could be explained by large distance 
modifications of our present understanding of gravity and inertia, without postulating non- 
baryonic forms of matter such as dark matter and dark energy [11], [12]. One of the most 
promising possibilities is represented, in the context of the brane-world scenario, by the 
Dvali-Garbadze-Porrati (DGP) model \T\. In the DGP model the gravitational interaction 
is described by five- dimensional Einstein gravity and the usual four-dimensional (4D) 
Newtonian potential is generated by a 4D scalar curvature term in the world-volume 
brane action. 

The idea of modifying Newton's law and/or Einstein gravity at large distances is not 
new. Since galactic rotational curves were found to be inconsistent with the distribution 
of luminous matter, such alternative theories of gravity have been proposed and pushed 
forward [21 [H, [6]. On the other hand, from the experimental side, there is plenty of room 
for such infrared modifications of gravity; our experimental knowledge of gravity conforms 
with general relativity but is limited to distances between say 10~^cm — lOOMPc. 
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A quite generic prediction of infrared-modified theories of gravity is the presence of 
large-distance logarithmic corrections to the Newtonian potential. In the DGP model, 
a ln(r/L) term, where L is a cross-over length, appears as correction of the Newtonian 
potential [U [T3]. In models that modify Newtonian dynamics to account for galactic 
rotational curves the logarithmic corrections dominates for small accelerations, of order 
Oo ~ 10~^^ cm~^ (in natural units) [31 H]. 

The search for Einstein-like theories of gravity that allow for a non-Newtonian weak- 
field limit with large-distance logarithmic correction is therefore of great interest. Such 
models can be considered as an effective descriptions of more general theories of grav- 
ity, such as those emerging in the brane-world scenario, which modify the gravitational 
interaction at large distance. 

It is not easy to build metric, Einstein-like, theories of gravity that allow for a non- 
Newtonian weak-field limit. Typically one has to introduce additional fields (scalars or 
vectors) and add to the Einstein-Hilbert action additional terms. Alternatively, one may 
consider higher-curvature terms in the gravitational action, e.g. in the form of f{R) 
theories. In both cases the theory is highly constrained by several phenomenological 
requirements: existence of a Newtonian potential term at short distances, absence of a 
fifth force, correct predictions for light bending and so on (see e.g. Ref. |i6j). 



A simple theory of Einstein gravity coupled to a scalar field with an exponential po- 
tential that allows for a ln(r/L) weak field limit has been proposed in Ref. [Hj. The 
proposed gravity theory requires a particular form of both the scalar field potential and 
the coupling function between scalar field and matter. 

In this paper we show that the existence of a ln(r/L) weak-field limit is not restricted 
to the particular model of Ref. [I^, but is a quite generic feature of a broad class of 
models of gravity coupled with a scalar field. We require that the model satisfies the 
following conditions: 

a) Existence of a consistent non-relativistic, static, weak-field limit of the gravitational 
field equations. 

b) At large distance (r ^ L) the gravitational potential is dominated by the ln(r/L) term. 

c) Presence of a sub-leading Newtonian term —GM/r. 

d) The short distance deviations from Newtonian behavior are at most of order r~'^ with 
7 > 2. 

We show that virtually all models of Einstein gravity coupled with a scalar field satisfies 
the above conditions if an appropriate coupling function between the scalar field and 
matter is introduced. As an illustration of our procedure we discuss in detail the case of a 
scalar field with a power-law potential. The models can be also rewritten, in the Jordan 
frame, as a scalar-tensor theory of gravity and also reformulated as a theory of gravity 
with an f{R) Lagrangian. 

The structure of the paper is as follows. In section 2 we present our model. In section 3 
we discuss the weak-field limit. In section 4 we derive the general solution. In section 5 we 
illustrate our general procedure for the particular case of a scalar field with a power-law 
potential. The relevance of our models for describing the infrared modifications of the 
gravitational interaction is discussed in section 6. In section 7 we formulate our model in 
the Jordan frame and we stress its equivalence with a f{R) theory of gravity. Finally, in 
section 8 we state our conclusions. 
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2 The model 



Let us consider a system of two point-particles of mass M (the source) and m (test 
particle), with M ^ m interacting with the gravitational field g^i, and a scalar field 
(p. The gravitational field is described by the Einstein-Hilbert action. The scalar field 
is self-interacting, with potential V{ip,L), and its interaction with the point-particles is 
characterized by a coupling function e^^'^'^\ Both the potential and the coupling function 
depend not only on ip but also on a real parameter L with dimensions of a length (we use 
natural units). For a generic, V{ip) ^ Ay^"^, potential the presence of such a parameter 
is necessary for purely dimensional reasons, because the scalar field has dimension of 
a (length)"^. We will show later in detail that this parameter has the physical meaning 
of the length-scale, above which the non-Newtonian behavior of our gravitational model 
become relevant. 

The classical dynamics of the system is described by the Einstein-like action 



5" = y d'^x^/^ 



1 



16nG 



R - d^ipd^ip - V 



j ds(^M) + m J ds 



(m) , 
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where ds[M,m) are respectively the line elements of the source and test particle ds(M,m) = 
dt[{dx1^j^ ^•^/dt){dx'^j^^f ,^^/dt)g^uY^'^ '^^ using for the metric a signature (—1, 1, 1, 1). 

In the limit M ^ m the contribution of the test particle to the stress-energy tensor 
and to the equation for the scalar field can be neglected. Moreover, we are only interested 
in the motion of the test particle. The resulting equation of motion stemming from the 
action ([1]) are, 



SnG 



d dX(^Yn)^ dX(^Yn) 

' per 



MH'e" I dr 



d^ix" - x" 



(M) 



V)) 



(2) 



dr^ 



dr dr 



0, 



where the prime denotes derivation with respect to (p> and T^^-*, T^^^ are the stress- 
energy tensors for the scalar field and for the source of mass M, respectively, and F is a 
(y9-dependent connection: 



p,u 

rp(M) 
fJ,U 



2d^(pd^ip 



M / drUnUy 



(3) 



In the previous equations is the four-velocity of the source. 



3 The weak-field limit 



We are interested in the usual weak-field, non-relativistic, static limit of the field equations 
([2]). We expand the metric near a flat background, g^y = rj^^ + h^^ with h^^ <^ 1 and we 
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consider field configurations depending only on the spatial coordinates x', i = 1,2,3. In 
the non-relativistic limit, when the velocity of the particles are <^ 1 and |Tjj| <^ |Too|, the 
field equations ([2]) give. 



AttG 



H 



12^ 



(V^)^ + V\+ H"{VvY + l^H'V + (^4vrG + — ) -oo 
1 

2' 



(M) 



(4) 



(fx 



(rn) 



-V<l>, 



where $ 



-hoo/2 + H, f^^'^ 



M6Hx 



X 



(M)j 



The bar indicates three-dimensional 



vectorial quantities, and the differential operators are calculated with respect to the three- 
dimensional Euclidean metric. 

From the third Equation in (jlj) it is evident that in the weak-field limit the field $ 
represents the potential that determines the force acting on the test particle. The usual 
weak-field Newtonian limit can be trivially recovered setting in Eqs. ^ H = V = Q 
and picking the </9 = solution for the scalar field equation. In the following we will 
only consider spherical symmetric solution to the Eqs. (jl]). They can be found placing 
the source particle of mass M at the origin of the coordinate system and using spherical 
coordinates ir,6,uj). 



4 The general solution 

Our goal is to find a general class of models described by the action ([T]) and parametrized 
by a particular, albeit general, from of the functions V^(v9, L), H{ip,L), whose solutions 
satisfy the following conditions: 

a) Existence of a consistent non-relativistic, static, weak field limit given by Eq. (jlj). 

b) At large distance (r ^ L) the solutions for $ are dominated by a ln(r/L) term . 

c) They must allow for a sub-leading Newtonian term —GM/r. 

d) The short distance deviations of $ from the Newtonian behavior are at most of order 
r~'^ with 7 > 2. 

Obviously, these conditions constraint the form of the scalar potential V and on the 
coupling function H . Consistency of the non-relativistic approximation requires that the 
stress-energy tensor satisfies the following condition: \Tij\ <C |Too|. In the case of the 
source of the gravitational field the previous conditions are satisfied if one considers a 
point particle moving with speed much lesser then the speed of light. On the other hand 
for the stress-energy tensor of the scalar field in Eq. ([3]) we have: 

Tl^^ = {drvf + V, T(^)'^ = -4^-2(9.^)2. (5) 

It follows that the conditions for the validity of the non relativistic approximation are 
consistent with the weak-field equations (jl]) only if V and {dj-ipY do not affect the New- 
tonian term in the potential This together with condition d) above imply the large r, 
r ^ L behavior: 

V{^) ~ (drif)^ ~ 0{r-^), (3>4. (6) 
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Condition c) requires the presence of the term —MG/r in the solution for $. In view 
of the first equation in (jl]) this imphes [AtiG + H''^)e^\r=o = 47rG. The simplest solution 
to this equation is given by 

H(r = 0) = H'{r = 0) = 0. (7) 

Our next and last task is to fulfil condition b). Using ([7]) and the equations of motion for 
the scalar field (p, the first equation in (jlj) becomes 

= V'ii + 47rGTo^'\ (8) 

In order to have a solution for the gravitational potential $ satisfying at the same time 
condition d), we must have for r ^ L 

W'H = ^ + 0{r~'), (9) 

where C is an arbitrary constant. Obviously Eq. determines the dependence of H on 
r. In order to find the dependence of the coupling function H on cp, we need to solve first 
the equation of motion for (/^(r). The r ^ L asymptotic form of the solution of Eq. (Q is 

H = C\ny + Ci + Oir-^), (10) 

where Ci is an arbitrary constant. Notice that in principle we can use Eq. ([8]) to describe 
an arbitrary deviation (not necessarily logarithmic) of the gravitational potential from its 
Newtonian behavior. Obviously in this case Eq. (Q as to be modified accordingly. 

We have now characterized completely the general model we are looking for. Every 
theory of gravity coupled with a scalar described by the action ([T]), will allow in the 
non-relativistic, static, weak-field limit for a solution 

r GM 

$ = Cln--^ + 0(r-2), (11) 
L r 

if the coupling function satisfies Eqs. fl71). f|T0l) . For instance a simple solution of Eqs. 

(IZD^dinD is 

^^ = fln(l + ^). (12) 

The form of the potential V{(p) is only constrained by Eq. ([6]), which is necessary 
for the consistency of the non-relativistic, weak-field approximation. An other physical 
condition that has to be imposed on the form of the potential is stability, i.e. V{ip) must 
be bounded from below and eventually it must allow for a local minimum. Excitations 
near this minimum must have enough heavy mass to be compatible with particle physics 
phenomenology. 



5 Scalar field with a power-law potential 

As a simple example and illustration of our general approach we consider the case of a 
scalar field with a power-law potential 

V = NL-\ipLY^-^/''\ (13) 
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where and a are arbitrary dimensionless parameters. The weak-field equations (jl]) for 
are readily solved to give 



where D = {[N{a — l)]/[a^(a + 1)] . Using solution (HM one easily find that the con- 
sistency conditions (Q are satisfied if a < —1. 

The coupling function H is constrained by Eqs. (17|) . (fTO!) . A simple solutions is given 
by Eq. (fT2|) . which in view of Eq. (HM becomes: 



For a < — 1 the potential V{!f) is bounded from below in the range < (f < oo, corre- 
sponding to the physical range < r < oo of the radial coordinate r. Thus the system 
is stable and in principle one can also deform the potential in such a way that V{ip) has 
local minimum at = ipo, with V"{(fQ) = m^. The parameters in the potential have to 
be chosen in such a way that the mass of the local excitation is heavy enough to be 
compatible with the experimental results about the absence of long range scalar forces. 

6 Large-distance modifications of gravity 

The gravity model ([T]) can be used both as an effective description of infrared modification 
of gravity, e.g. in the brane world scenario, or also as an Einstein-like description of 
Newtonian modified dynamics. 

In the first case the scalar field ip could represent an effective 4D parametrization of 
those effects that modifies gravity at large distances ( for instance the embedding of our 
3-brane in a higher dimensional world). Apart from the potential V{(p), which does not 
determine the infrared modifications of gravity, the model is parametrized by a coupling 
function H which depends on a length scale L and a coupling constant C. H is zero 
together with its first derivatives at r = and stays almost constant (and vanishing) for 
r ^ L. In this regime the gravitational potential is given by the Newtonian expression 
and the theory reduces to General relativity. For r ^ L, if ~ ln(r/L), the gravitational 
potential is dominated by a In r term and the theory deviates from general relativity, still 
giving an effective metric description of the gravitational interaction. 

In the second case, one can use our model to explain the rotation curves of the galaxies 
without postulating the presence of dark matter. To this end one introduces a constant 
acceleration oq ~ 10~^^cm such that for a ^ one recovers standard Newtonian dy- 
namics. In fact, in this case yGM /r ^ cq and the Newtonian term in the potential (ITTi) 
dominates over the logarithmic term. On the other hand for a ~ Oq the leading loga- 
rithmic term dominates the gravitational potential ffTTl) and the rotation curves of the 
galaxies can be explain using the gravitational potential (fTTl) and identifying the constant 
C as 



where M is the mass of the galaxy and a is an universal constant, which can be determined 
using the Tully-Fisher law [15], [H] . 




(14) 




2 " 



(15) 




(16) 
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7 Jordan frame and equivalence with f{R) theories 

The gravitational model described by the action ([T]) is equivalent to a scalar-tensor theory 
of gravity. It is a scalar-tensor theory of gravity written in the Einstein frame. The 
transition from the Einstein to the Jordan frame is given by the Weyl transformation, 

g^Lu = e~'^"g^^, (17) 

where g^j^y is the metric in the Jordan frame. In this frame the action ([1]) becomes 



S = Id xJ—ge 



-2H 



^ Ma) + (^H' - 1 ) d^ifid'^ip - e-2^r 



+ 



IGttG \87rG 

- [Ml dsM + m f dsrr^ . (18) 



In the Einstein frame theories described by the action ([T]) lead to a breakdown of the 
equivalence principle at galactic scales [13] . Conversely in the Jordan frame the particles 
do not couple to the scalar field. The geodesic equation depends only on the metric g^,^. 
On the other hand in the Jordan frame the theory has the usual weaknesses of scalar- 
tensors theories of gravity: Newton constant depends on the coordinates and light bending 
is not correctly reproduced (see e.g. Ref. [H]). 

One can also consider the weak-field limit of the field equations stemming from the 
action f|T8l) along the lines described in the previous sections. Because now the geodesic 
equation for the test particle depends only on the metric g^^, we set g^^i; = rj^y + h^y and 
$ = —hQQ/2. In the non-relativistic, weak-field limit we can approximate e^^^ = 1 it 2H 
and the equation of motion for $, ip and the geodesic equation coincide with (j4]) upon 
use of Eqs. ([n]),(IZl)- Hence the solutions (ITU]) . flTT]) . f[T^ found in Sect. 4 are also solution 
of the theory in the Jordan frame. 

It is well known that gravity models whose Lagrangians are given by functions of the 
scalar curvature f{R) are on-shell equivalent to Einstein gravity coupled to a scalar field 
(see e.g. Ref. [6]). General models may also involve a an arbitrary coupling function e^^'^^ 
[H] and be therefore equivalent to the model described by the action ([T]). The potential 
V{!f) for the scalar field is determined by the the function f{R)- Our characterization of 
gravity models with non-Newtonian weak-field limit can be therefore easily extended to 
the case of gravity theories with f{R) Lagrangians. 

The existence of non- Newtonian solutions seems a quite generic feature of f{R) theories. 
In fact the potential V is determined entirely by the function / and the requirement for 
the existence of non-Newtonian solutions constrain very weakly the form of V. Also in 
this case the existence of the non-Newtonian weak-field limit is related to the presence of 
a non trivial coupling function . 



8 Conclusion 

In this paper we have show that the existence of a non-Newtonian weak-field limit in 
which the gravitational potential has a large distance ln(r/L) leading term is a quite 
generic feature of a broad class of models of Einstein gravity coupled with a scalar field 
and of theories of gravity with an f{R) lagrangian, when a scalar field- dependent matter 
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coupling function is introduced. This class of models has been selected by enforcing rather 
general physical conditions on the form of the non-relativistic gravitational potential. 

The gravity models proposed and investigated in this paper may be very useful to give 
an effective description for large-distance deviations of the gravitational interaction from 
its GR behavior, e.g in brane-world scenarios. Another interesting field of application is 
their use to describe modifications of Newtonian dynamics at small accelerations, e.g. for 
explaining rotation curves of galaxies. 

The most sensible and critical point of our approach is the physical origin of the 
coupling function 6^^"^^ Its presence is quite natural in scalar-tensor theories of gravity, 
e.g. Brans-Dicke theory and related to the choice of what may be called the physical Weyl 
frame [161 [13 • But its physical interpretation is far from being clear. One possibility is 
that it gives a local effective description of non-local effects. For instance, in brane-world 
scenarios it may be considered as a local parametrization of nonlocal effects generated by 
the embedding of the 3D brane in the higher dimensional spacetime. 
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